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I. INTRODUCTION 

The BaBAR detector [1], designed to study time- 
dependent CP-violation in B decays, is very well suited 
to study hadronic final state production. At an e~^e~ col- 
lider, lower center-of-mass (cm.) energies can be studied 
using initial state radiation (ISR). The possibility of ex- 
ploiting such processes was first outlined in Ref. [2] and 
discussed in Refs. [3-5]. 

The cross section a^f to radiate a photon of energy Ej 
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and subsequently produce a definite hadronic final state 
/ is related to the corresponding e'^e~ production cross 
section cr/(s) by: 



du^f{s, x) 



dx 



W{s,x) •cr/[s(l-x)] 



(1) 



where 



s, ^/s is the nominal e+e cm. en- 



ergy and Ec.m. = \^ s(l — x) is the effective cm. energy 
at which the final state / is produced. This quantity 
can be determined, for instance, by measuring the in- 
variant mass of the hadronic final state. The radiator 
function W{s,x), which describes the probability den- 
sity for photon emission, can be evaluated with better 
than 1% accuracy (see for example Ref. [6]). The radi- 
ator function falls rapidly with increasing Ej, but has a 
long tail, which produces sizeable event rates at very low 
Ec.-ax.- In the present study we require the ISR photon 
to be detected in the BaBM. electromagnetic calorimeter. 
The angular distribution of the ISR photon peaks along 
the beam direction, but 10 — 15% [5] of the photons are 
emitted within the detector acceptance. 



II. THE BABMl DETECTOR AND DATASET 

The data used in this analysis, collected with the 
BaBAR detector at the PEP-II storage ring, correspond to 
an integrated luminosity of 232 fb~^. Charged-particle 
momenta are measured in a tracking system consisting 
of a five-layer double-sided silicon vertex tracker (SVT) 
and a 40-layer central drift chamber (DCH), immersed 
in a 1.5 T axial magnetic field. An internally reflecting 
ring-imaging Cherenkov detector (DIRC) with quartz bar 
radiators provides charged-particle identification. A Csl 
electromagnetic calorimeter (EMC) is used to detect and 
identify photons and electrons. Muons are identified in 
the instrumented magnetic flux return system (IFR). 

Electron candidates are selected using the ratio of the 
shower energy deposited in the EMC to the measured mo- 
mentum, the shower shape, the specific ionization in the 
DCH, and the Cherenkov angle measured by the DIRC. 
Muons are identified by the depth of penetration into 
the IFR, the IFR cluster geometry, and the energy de- 
posited in the EMC. Charged pion and kaon candidates 
are selected using a likelihood function based on the spe- 
cific ionization in the DCH and SVT, and the Cherenkov 
angle measured in the DIRC. Photon candidates arc de- 
fined as clusters in the EMC that have a shape consistent 
with an electromagnetic shower, without an associated 
charged track. 

In order to study the detector acceptance and effi- 
ciency, a special package of simulation programs for ra- 
diative processes has been developed. Algorithms for 
generating hadronic final states via ISR are derived from 
Ref. [7]. Multiple soft-photon emission from initial- 
state charged particles is implemented by means of the 
structure- function technique [8, 9], while extra photon 



radiation from final-state particles is simulated via the 
PHOTOS package [10]. 

The e+e~ — > K'^gK^ir^^, A'+if"7r°7, 4>ri^ and 07r°7 
reactions are simulated with the hadronic final states gen- 
erated according to their phase space distributions. Both 
> P1P2P3 and e+e~ VP [7], where 



processes e^e 



and V arc pseudoscalar and vector mesons, respec- 
tively, are used in the simulation. To evaluate the back- 
grounds, a number of different ISR-produced final states, 
such as 47r, KKtttt, and 07r°7r°, are simulated as well. 

The quark-antiquark continuum production and 
hadronization (qq background in the following, where 
q = u,d,s,c) and rr events are simulated using JET- 
SET [11] and KORALB [12], respectively. 

The generated events undergo a detector simulation 
based on GEANT4 [13] and are analyzed as experimental 
data. 



III. EVENT SELECTION AND KINEMATICS 

A preliminary set of selection criteria are applied to 
the data regardless of the final state under study. We 
consider only charged tracks emitted at a polar angle be- 
tween 0.45 and 2.40 radians and that extrapolate within 
0.25 cm of the beam axis and within 3 cm of the nominal 
collision point along the axis. We retain photon candi- 
dates with energy above 30 MeV and emitted at a polar 
angle between 0.35 and 2.40 radians. The highest-energy 
photon, assumed to be the ISR photon, is required to 
have an energy greater than 3 GeV. 

We then apply channel specific selection criteria and 
subject each candidate event to a one-constraint (IC) 
kinematic fit, both to select the signal and to estimate 
the background contamination. The fit uses as input 
the measured quantities (momenta and angles) of the 
selected charged tracks and photons, together with the 
associated error matrices. However, because of the excel- 
lent resolution of the DCH, the three-momentum vector 
of the ISR photon is better determined through the mo- 
mentum conservation than through the measurement in 
the EMC, thus reducing the number of constraints to one. 
For final states containing a tt^ or an rj, we perform the 
kinematic fit for every possible pair of detected photons, 
without applying a mass constraint, and retain the pair 
producing the lowest x^. The fitting procedure predicts 
the ISR photon direction, which is then required to be 
consistent with the observed ISR photon within 15 mrad. 
In order to reduce multi-photon events, we also require, 
that no more than 400 MeV of additional neutral energy, 
besides the ISR and 7r° (77) decay photons, be detected 
in the EMC. This cut has no effect on signal events. 

We select as 7r° (77) candidates 77 pairs whose masses 
lie in the range 0.110 - 0.160 (0.520 - 0.580) GeV/c^. 
In calculating the invariant masses of these pairs, we 
use the photon momenta resulting from the IC fits de- 
scribed above. This improves the tt*^, rj mass resolution 
by w 20%. To get rid of the contamination due to addi- 



7 



tional soft photons produced by machine background or 
by interactions in the detector material, we require both 
photons used for tt" or rj candidates to have an energy 
greater than 60 MeV and at least one cluster containing 
more than 100 MeV. 

As Kg candidates, we select pairs of charged tracks 
forming a vertex with a good least 2 mm away 

from the event primary vertex, and whose tt+tt" invari- 
ant mass lies within 15 MeV/c^ of the nominal Kg mass. 
We also require that d^o, the angle between the Kg 
momentum vector and the line connecting the Kg ver- 
tex to the primary vertex position, satisfy the condition 
cos(6'/fo) > 0.99. For the purposes of the IC fit described 
earlier, the Kg candidate is treated as a single particle 
and its covariance matrix is used as the corresponding 
track covariance matrix. 

Charged kaon and pion candidates are selected using 
likelihood ratios for the hypotheses e, ^, tt, K, and p 
based upon specific ionization measurements in the SVT 
and DCH, and Cherenkov light observed in the DIRC. 
The algorithm used here for K^ identification is over 
80% efficient in the momentum range of interest, with tt, 
ji, p, e ^ K misidentification rates well below 2%. 




10 20 30 40 

FIG. 1: The x'^ distributions for data (points), signal MC 
(open histogram) normalized to the first bin of the data, back- 
ground contribution due to qq (hatched histogram) and dif- 
ferent ISR-produced final states (gray histogram) . 



IV. THE KlK^-K^ FINAL STATE 
A. Event selection 

We require exactly four charged tracks with zero net 
charge, two of them coming from the primary vertex 
of the event and identified respectively as a kaon and a 
pion, and the other two forming a _ft'". 

Figure 1 shows the distribution of the selected 



data sample, together with three different MC sim- 
ulations: signal events, background contribution 
due to e^e~ — s- qq, and ISR-produced processes 
K'^K~TT^'K~ + 2(7r+7r~) that might mimic this final 
state. A total of 3860 events are selected within the x^ 
signal region, i.e. x^ < 20. 

The background due to mis-reconstructed K^ is small, as 
can be inferred from the K^ mass distribution shown in 
Fig. 2, and it is almost exclusively given by the reactions 
e+e^ — > ^K'^ K~'K'^TT~ and e^e~ — *■ 77r"'"7r~7r"'"7r~. 
These processes differ from the signal only by one K — tt 
exchange, and produce a x^ distribution peaking at low 
values. The direct simulation of these two sources of 
background yields to respectively 44 ± 15 and 10 ± 5 
events in the x^ signal region, where the errors are taken 
from the uncertainties of the corresponding production 
cross sections [14]. 
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FIG. 2: Reconstructed A'g mass for data (points) and signal 
MC simulation (histogram). 

The cross sections for the background processes in 
continuum qq production are poorly known; this type of 
background is almost exclusively produced by processes 
like e+e" — *■ K'^K^tt^tt'^ with the 7r° decaying into two 
very energetically asymmetric photons and/or photons 
merged in a single EMC cluster that fake the ISR pho- 
ton. This class of events, as a consequence, produces a 
sharp peak at the tt'^ invariant mass when the fake ISR 
photon is combined with another photon in the event. 

Figure 3(a) depicts the invariant mass combination of 
the candidate ISR photon with any other photon in the 
event. A sharp peak at the 7r° mass is evident for the qq 
MC sample, while the signal and the data distributions 
are both essentially flat in that region. We perform a 
fit to the data distribution, combining with weights the 
distributions of qq background and signal. The difference 
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M^(GeV/c') 



M^(GeV/c-) 



As a cross-check, we show in Fig. 3(b) the invariant 
mass distribution of any pair of detected photons in the 
event, not including the ISR photon. The very good 
agreement between the data and the simulated signal 
sample and the absence of a peak at the tt" mass, con- 
firm that background contribution from final states with 
extra 7r°'s is negligible. 

The K^K^TT^ invariant mass distribution for the final 
sample is shown in Fig. 4, together with the estimated qq 
and ISR backgrounds. The final signal yield after back- 
ground subtraction consists of 3714 events. 



FIG. 3: (a) invariant mass of the ISR photon and any other 
photon in tlie event. Distributions for data (points with er- 
ror bars), signal MC (open histogram) and qq MC (shaded 
histogram) samples are shown, (b) 77 mass distribution for 
additional photons in data (points) compared with signal MC 
simulation. 



between the shapes of the signal and qq distributions al- 
lows to obtain unambiguously the weights from the fit. 
We find a weight for the qq of about 1/9, which translates 
into 72 ± 42 qq background events. 

Any remaining background due to other ISR processes 
should results in a relatively flat distribution, as a con- 
sequence of the mis-reconstructed final state. We check 
the consistency of the distribution by normalizing the 
data and the simulated K'^K^n^ sample to the first bin 
of the data distribution (see Fig. 1). after subtracting 
the already estimated background components. The dif- 
ference between data and normalized signal distributions 
gives an estimate of additional background events con- 
sistent with zero. 




M(K„K7i:) (GeV/cl 



'^TT^ sys- 



FIG. 4: The invariant mass distribution of the KgK 
tem (open histogram). The hatched and gray histograms 
show the qq and ISR backgrounds, respectively. 
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FIG. 5: Detection and reconstruction efficiency as a function 
of the cm. energy for the K^K'^n^ final state. The solid 
line is the result of a linear fit to the data. 



B. 



Cross section for e~'~e 



The cross section for the e+e" K2,K^tt^ process as 
a function of the effective cm. energy, i^c.m.j is evaluated 
from: 



[Ec.-m.) 



(2) 



where Ecm. is the invariant mass of the reconstructed 
K'^K^tt^ system, dNj^oj^^^ is the number of selected 
K^K^TT^ events after background subtraction in the in- 
terval dEc.m., and e{Ec.m.) is the corresponding detection 
efficiency obtained from the signal MC simulation. The 
differential luminosity, dC{Ec.m.) , in each interval dEc.m., 
is evaluated from ISR fifx'j events with the photon in 
the same fiducial range used for simulation; the proce- 
dure is described in Refs. [15, 16]. From data-simulation 
comparison, we conservatively estimate a systematic un- 
certainty of 3.0% on dC This error also includes the 



uncertainties on radiative corrections, estimated from a 
MC simulation to be sa 1%. 

The efficiency, that accounts also for the K'^ tt+tt^ 
decay rate, has been corrected for differences between 
data and simulation: the angular distributions are differ- 
ent from the simulated distributions because the latter, 
based on a phase space model, do not account for the 
presence of intermediate resonant states. However, the 
angular acceptance is essentially uniform for ISR events, 
so the uncertainty in estimated efficiencies due to impre- 
cise modeling of the process dynamics is, at most, 3.0%. 
We study several control samples to determine any differ- 
ence on track reconstruction efficiency between data and 
simulation. We find that the MC efficiency is overesti- 
mated by 0.8% per charged track, with an uncertainty of 
±0.5%. Analogously, we determine a correction of 2.6% 
on the A'g reconstruction efficiency, with an associated 
systematic uncertainty of 1.1%. The uncertainties on qq 
and ISR background estimates lead to systematic errors 
of 1.1% and 0.5%, respectively. 

The main systematic uncertainties assigned to the 
measurement are summarized in Table I. The corrections 
applied to the measured cross sections sum to -1-4.2%, 
while the estimated total systematic error is ±5.1%. 
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FIG. 7: Comparison of the e e ix^i 
measured in BABAR with previous experiments, in the mass 
range 1.3 < Sc.m. < 2.4 GeV. 
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FIG. 6: The e^e K^K'^n^ cross section. Inset: ex- 
panded view in the mass range 3.2 < iJc.m. < 4.6 GeV. 

For a typical ISR reaction, the selection efficiency is 
expected to have a smooth behavior as a function of 
Ec.m.- Therefore, we fit the efficiency values obtained 
from MC as shown in Fig. 5, and use the fit result, to- 
gether with the estimated corrections, to obtain the cross 
section plotted in Fig. 6. The cross section has a maxi- 
mum of about 3.8 nb around 1.65 GeV, and a long tail. 
The only clear structure observed at higher energies is 



TABLE I: Summary of corrections and systematic errors ap- 
plied to the e+e^ ^ KsK^tt^ cross section measurement. 



Source 



Correction Systematic error 



ISR luminosity 




3.0% 


qq background 




1.1% 


ISR background 




0.5% 


Track reconstruction efficiency 


+1.6% 


1.0% 


PID efficiency 




2.0% 


Ks reconstruction efficiency 


+2.6% 


1.1% 


B{K° ^Tv+n-) [17] 




0.7% 


MC model 




3.0% 


Total 


+4.2% 


5.1% 



the J/ip . In the inset of the same figure the cross section 
at energies higher than the J/tJj is shown: there is no ev- 
idence of the recently discovered y(4260) resonance [18], 



E^,^(GeV) we obtain an upper limit for F 



y(4260) „F(4260) 



of 0.5 eV 

at 90% C.L.. An excess of events is observed at « 4.2 
GeV: there are 15 events in a 40 MeV energy window, to 
be compared to an expectation of 5 events, corresponding 
to a probability of 2.2x10""' to be a statistical fluctua- 
tion of the cross section. 

The values of the measured cross section are reported 
in Table II with the corresponding statistical errors. Fig- 
ure 7 shows a compilation of the measured cross section 
as obtained by previous direct e'^e~ experiments (data 
are available only up to 2.4 GeV [19]). The data show 
reasonable consistency; our measurements, however, have 
much smaller statistical and systematic errors. 
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TABLE II: Measurement of the e+e" KIK^ TT^ cross section as a function of Ec.m.- Errors are statistical only. 
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V. THE K^R-'K^ FINAL STATE 

We select K^K^-k^ candidate events by requiring two 
charged tracks identified as kaons, and two photons form- 
ing a tt" candidate, in addition to the ISR photon. In 
order to avoid any possible contamination from events 
with intermediate — s- K^K~ production, we require 
^KK > 1.045 GeV/c^, where Mkk is the invariant mass 
of the K^K^ system. As previously described, we per- 
form the kinematic fit for every possible pair of selected 
photons and we retain for each event the 77 combination 
giving the lowest . Figure 8 shows the distribution 
of the selected data sample, and the normalized distribu- 
tions of signal MC simulation and different background 
sources {qq and e"''e~ lo^ 7r+7r^7r*'7). The selection 
procedure results in 1524 events with a < 20. 

Applying a procedure analogous to the one described 
in the previous section for the K^K^i:^ analysis, we de- 
termine 62 ± 35 qq and 84 ± 27 ISR-background events 
from different ISR processes. The mass spectrum of the 
selected sample is shown in Fig. 9, together with the es- 
timated qq and ISR backgrounds. After background sub- 



traction, 1378 signal events are retained. 

As in the K^K^tt^ channel, we measure the e^e'^ 
K~^K~7t'^ cross section as a function of the effective cm. 
energy, in 20 MeV energy bins. The reconstruction and 
selection efficiency for each energy interval is shown in 
Fig. 10; while the corresponding cross section values are 
shown in Fig. 11 and reported in Table III (variable width 
bin size). Figure 12 shows a comparison with the DM2 
measurement [19]. 

The cross section values shown in Fig. 11 include only 
statistical errors. The inset of the same figure is an en- 
larged view of the 3.3 — 4.6 GeV energy region, there 
is no evidence of the y(4260), and we estimate an up- 
per limit for rri''^^"'S^^^^^°^o of 0.6 eV at 90% C.L.. 
The systematic errors are summarized in Table IV, along 
with the corrections applied to the cross section measure- 
ments [15, 16]. The correction and systematic error asso- 
ciated with the MC /data difference in the tt" reconstruc- 
tion efficiency have been studied using the ISR-produced 
high statistics channel wtt" 7r"'-7r~7r°7r° [16]. It has 
been found that the MC efficiency is 3.0 ± 3.0% higher 
than the data. The corrections applied to the measured 
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40 

X'(K"K-7U°) 




M(K+K"7U°) (GeV/c^) 



FIG. 8: The x distributions for data (points), signal MC 
simulation (open histogram) normalized to the first bin of 
the data, background contributions due to qq (hatched his- 
togram) and to e^e" luj 7r^7r~7r''7 MC simulation (gray 
histogram) . 



FIG. 9: The invariant mass distribution of the K'^K~-k'^ sys- 
tem (open histogram). The hatched and gray histograms 
show the qq and ISR backgrounds, respectively. 



The Dalitz method 



cross sections sum to -1-4.6%, 
systematic error is ±6.5%. 



while the estimated total 



VI. THE DALITZ PLOT ANALYSIS FOR 

° CHANNELS 



/fsA:±7r=F AND K+K n 



Figures 13 and 14 show the Dalitz plot at cm. en- 
ergies below 3.1 GeV for the K+K~Tr° and K^K'^tt^ 
final states, respectively. We observe that the in- 
termediate states KK* (892) and ^^^(1430), carrying 
the allowed combinations of quantum numbers, give 
the main contributions to the KKn production, while 
KKl{14:10) contributes only weakly due to the small 
value of r{K^{U10) ^ Kn). 

The final state K~^K~tt'^ can be produced through 
the intermediate states K*'^{892)K'f' or K*^{U30)K'^, 
and we expect a Dalitz plot with a symmetric popula- 
tion density with respect to the exchange K~^Tr^{K*'^) ^ 
K'-K^iK*-), see Fig. 13. 

The final state K^K^tt^ is obtained via the decays 
oiK*^{892)K'F and K*°{892)K°, or K;^{U30)K'F and 
Ki°{U30)K^, I.e., both charged K*^{K°n^) and neu- 
tral K*°{K^tt'^) can be produced. The population den- 
sity of the Dalitz plot is expected to be asymmetric in 
this case. This effect is clearly seen in Fig. 14. 

A detailed explanation of the analysis techniques and 
the results obtained is given in the following subsection. 



The Dalitz plot population has been fit assuming KK* 
two-body processes as intermediate states and a contin- 
uum contribution coming from tails of resonances cou- 
pled with KK* , but lying below their production thresh- 
old [20] (see Sec. X). The density function is written in 
terms of 3- vector isospin amplitudes Ai_k* , that, in turn, 
are products of a Breit Wigner distribution {BW), the 
K* propagator, and angular distribution probabilities, 
described by Zemach tensors [21]. In this case, since the 
spectator particle is always a pseudoscalar meson (the 
kaon), the Zemach tensors are simply 3- vectors and the 
amplitudes are: 

Al,K'=Ci.K'fK' {Prr - PK , ^k) ■ BWk' [{p^ + PK?] + 

(3) 

^^iji'^K' (Pk-P^iPk) ■ BW-j^'[{pj^ + p^f], 

where p^ is the 4-momcntum of the generic final meson a 
and Tfc* is the Zemach tensor related to the K*. The two 
terms of Eq. 3 correspond to the two possible Kit and Ktt 
combinations in which the K* and K can exist. The 
coefficient , which accounts for the total isospin / and 
the charge of the intermediate kaon a, has the form 



C 



7, if* 



(-1) 



1)^ for charged K* 
for neutral K* 



(4) 



The Zemach tensors for K*{892) and the K^{U30) are 
shown in Table V. 
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TABLE III: Measurement of the e^e — > K vr cross section as a function of -Bc.m.- Errors are statistical only. 
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TABLE IV: Summary of corrections and systematic errors 
applied to the e^e~ — > K~n^ cross section measurement. 



Source Correction Systematic error 



ISR luminosity 




3.0% 


qq background 




2.5% 


ISR background 




2.0% 


Track reconstruction efficiency 


+1.6% 


1.0% 


PID efficiency 




2.0% 


tt" reconstruction efficiency 


+3.0% 


3.0% 


MC model 




3.0% 


Total 


+4.6% 


6.5% 


TABLE V: Values of the Zemach 


tensor Tk' 


[20] in the two 


cases under consideration. 






Meson K* L 


TK'ipi 


-P2,P3) 


K'{892) 1- P (Pi-P2)xp3 
^-5(1430) 2+ D [(pi - pa) • PaWi - P2) x pz 



The Dalitz density has been fit in each bin of i?c.m. with 



the function 

V[{p^+PK)^,{pj^ + p^f,E,.„,]= (5) 

/,/'=0,l 

K* , K*' = K* (892), (1430) 

where the Aj^k* are the known amphtudes of Eq. 3, 
which depend on the K* (Kir) invariant mass and the 
total energy i?c.m.j and the complex coefficients fi^K-, 
depending only on Ec.m., and proportional to the isospin 
component of the corresponding K*K channel, are the 
free parameters. The complex values of the fi^K* pa- 
rameters have been used to extract moduli and relative 
phase of the isospin components. The analysis described 
in the following applies only to the K'^K^tt^ Dalitz plot, 
since only in this case, with both charged and neutral 
K* contributing, we can separately extract the isoscalar 
and isovector components. Because the KK* (892) and 
if ^2 (1430) represent the dominant contribution to the 
KKtt production at cm. energies below and above 
2.0 GeV, respectively, these two regions have been con- 
sidered separately. 

Figure 15 shows, in the low-energy region, the 
K'^K^tt^ Dalitz plot and the three possible projections. 
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FIG. 10: Detection and reconstruction efficiency as a func- 
tion of the cm. energy for K'^K^n'^ final state. Tlie 
solid line is the result of a fit to the data with the function 



The dominance of the i^*(892) in each Ktt channel ap- 
pears evident. As a result, in the following we will 
consider the process e+e^ K*{^%2)K as the domi- 
nant one. The different shapes observed for the charged 
(Fig. 15(b)) and neutral i^*(892) (Fig. 15(c)), are due to 
interference effects between the isoscalar and isovector 
components (see Sec. XA). 

Figures 16 and 17 show the isoscalar and isovector 
cross sections; numerical values are reported in Tables VI 
and VII respectively. Figure 18 shows the relative phase 
between the two amplitudes. In addition to the errors 
reported in Table I we must consider the 15% system- 
atic uncertainty associated with the Dalitz plot-fit pro- 
cedure. This quantity has been evaluated by studying the 
effects of different parameterizations for the BW distri- 
butions which describe the intermediate K* resonances 
in the amplitudes of Eq. 3. The isoscalar contribution, 
see Fig. 16, is dominant with respect to the isovector, 
(Fig. 17) and clearly shows a resonant behavior. More- 
over, the isovector cross section seems to be compatible 
with a pure phase space shape. However, using the com- 
plete information, including the relative phase, shown in 
Fig. 18, we also find that the isovector component has a 
resonant structure (for details see Sec. XA). 

Figure 19 shows the Dalitz plot and the projections 
obtained in the cm. energy region above 2.0 GeV. In 
this case, in addition to the K K* {%92) ^ we have evidence 
that the ifi^2(1430) intermediate state also contributes. 
In particular, as theoretically predicted [22], we observe 
that while the neutral i^*(892) dominates the tt^ 



FIG. 11: The e'^e — > K tt" cross section. Inset: ex- 
panded view in the 3.2 < Ec.m. < 4.6 GeV mass range. 




2 2.25 

E_(GeV) 



FIG. 12: Comparison of the e^e K tt" cross section 

measured in BABAR with previous experiments. 



projection (Fig. 19(c)), the charged ii'|(1430) is the main 
contribution in the tt^ projection (Fig. 19(b)). This 
can be interpreted as a different spatial structure of these 
mesons: _ftr*(892) is a vector, while X2(1430) is a tensor. 

The Dalitz plot analysis in this energy range is per- 
formed with four isospin amplitudes, an isoscalar and 
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TABLE VI: Isoscalar component of the cross section as a function of Ec.m. for the process e^e — > K* {S92)K . Errors are 
statistical only. 
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TABLE VII: Isovector component of the cross section as a function of Ec.m. for the process e^e —> K* {892)K . Errors are 
statistical only. 



-Ec.m.( 


GeV) 




a 


(nb) 


Sc.m.(GeV) 




a 


(nb) 


Ec.m. 


(GeV) 


a 


(nb) 






(GeV) 




a 


(nb) 


1.40- 


■1.42 





.32 


± 0.67 


1.64-1.66 


1, 


.83 


± 0.37 


1.88 


-1.90 


0.59 


± 0. 


,26 


2.28 


-2.32 


0, 


.25 


± 0.09 


1.42- 


■1.44 





.68 


± 0.40 


1.66-1.68 


2, 
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± 0.50 


1.90 


-1.92 


0.26 


± 0. 
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-2.36 


0, 
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.40 


± 0.43 
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an isovector component for each intermediate two-body 
channel. Figures 20(a) and 20(b) show the isoscalar and 
isovector components of the A'2 (1430)/^ channel. The 
corresponding values are reported in Tables VIII and IX 
respectively. The structure observed at Ec.m. — 2.1 GeV 
in the isoscalar and Ec.m. — 2.25 GeV in the isovector 
components may be interpreted as pure threshold effects. 
Indeed, the transition form factors, obtained by taking 
the square root of the cross sections divided by the ap- 
propriate phase space (a pseudoscalar and a tensor meson 
in the final state), have a steep yet smooth behavior, as 
shown in Fig. 21. The simplest interpretation calls for 
transition form factors that collect the main contribu- 
tions from low-energy resonances lying below threshold. 
Further study and higher statistics would help to clarify 
this item. 



VII. THE J/ip REGION 

Clean J/^ signals are observed in almost all the fi- 
nal states under study, from which we extract the corre- 
sponding decay rates. Figure 22(a) shows the inclusive 



distribution in the J/tp mass region for the K^K^tt^ 
data sample. A fit with a Gaussian for the J/ip peak and 
a linear function for the continuum, gives a{Mj^oj^^) = 

6.5 McV/c^ and less than 1 MeV/c^ difference from the 
PDG [17] value for the J/ip mass. To account also for 
possible radiative or experimental effects we determine 
the yields under the J/ip peak by fitting with a double- 
Gaussian. As shown by the Dalitz plot analysis, the 
e+e^ KKt: reactions proceed mainly through for- 
mation of the intermediate KK* (892) state. We could 
therefore extract the J/tp K*{892)K decay rates, 
where the K*{892) is charged {K*^{892) ^ K^tt^) for 
the K+K^n" final state and neutral i^*°(892) ^ K^tt^ 
or charged (i^*±(892) ^ K°tt^) for the K^K'^tt^ final 
state. In order to isolate the different if* (892) contribu- 
tions, we require Ek > 1.4 GeV in the KKtt cm. frame 
for the kaon not coming from the if* (892) decay, as ex- 
pected for a two-body J/tp decay. A check of the Dalitz 
plot provides an estimate of the residual contamination 
from wrong Kn combinations, which is found to be less 
than 3.0%. From the fit to the mass distributions ob- 
tained for the different A"* (892) A' final states, shown in 
Fig. 22(b-d), we derive the number of events reported in 
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FIG. 13: The Dalitz plot distribution for the A'+A'^tt" final 
state. 
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FIG. 15: (a) Dalitz plot distributions for A^A^Tr^ events 
with an invariant mass m^oj^^ < 2.0 GeV/c^. (b) Ag tt* 

projection, with the broad charged A* (892) peak, (c) A^ tt^ 
projection, with the narrow neutral A* (892) peak, (d) Ag 
A^ projection. 
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FIG. 14: The Dalitz plot distribution for the K°sK^n^ final 
state. 



Table X, together with the corresponding J/ip branch- 
ing fractions multiplied by the J/ip —^e^e" width. The 
latter quantity has been obtained through 



f ee 



where B'^^^'^ is the branching fraction for J/i]j /, 

dC/dE = 65.6 nb^VMeV is the ISR luminosity inte- 
grated at the J/ip mass, e is the detection efficiency ob- 
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FIG. 16: Isoscalar component of the cross section for the pro- 
cess e+e~ A* (892) A. The vertical line indicates the sep- 
aration between the two regions used in the analysis. 



taincd from simulation, and C = {hc/2'K)'^ = 3.894 x 
10^^ nbMeV^ is a conversion constant. The system- 
atic errors have been estimated by considering, in ad- 
dition to the values reported in Tables I and IV, the 
statistical error on the reconstruction efficiency. As- 
suming the world average value of the J/ip e^e~ 



width, r, 



5.55 ± 0.14 keV [17], we extract the decay 
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FIG. 17: Isovector component of the cross section for the 
process e^e~ K* {S92)K . The vertical hne indicates the 
separation between the two regions used in the analysis. 
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FIG. 18; Relative phase between the isoscalar and isovector 
component of the e^e~ K*{892)K cross section. 
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FIG. 19: (a) Dalitz plot distributions for KsK^tt^ events 
with an invariant mass m^o^Tr > 2.0 GeV/c^. (b) tt* 
projection, with the broad charged K2{U30) peak, (c) 
TT^ projection, with the narrow neutral if* (892) peak, (d) 
Kg projection. 




rates shown in the last column of Table X. Here and 
in the following the first uncertainty reported is statis- 
tical and the second is systematic. From these results 
and the known /'ir*(892) decay rates, we can derive the 

J/V' ^ K*°{892)lf,K*+{892)K- branching fractions, 
summarized in Table XI. 



VIII. THE 07r" FINAL STATE 

The event selection for (/)7r° production closely follows 
the procedure described for A'+A'^tt" with non-resonant 
A'+AT^pairs. We tighten a few cuts in order to remove 
random photon combinations that could fake a 7r°. In 
particular we lowered to 300 MeV the cut on the addi- 
tional detected neutral energy, and we require that the 
direction of the observed ISR photon matches that pre- 



FIG. 20: Isoscalar (a) and isovector (b) components of the 
cross section for the process e^e~ — > A"2(1430)A. 



dieted by the one-constraint fit within 12 mrad. 

Figure 23 shows the K'^K" mass, Mkk, versus the 
77 mass, M^-,,. The plot is divided into nine equally 
sized regions, with the central one, which corresponds 
to 1.010 < Mkk < 1-030 GeV/c^ and 0.110 < M^^ < 
0.160 GeV/c^, being the signal region. We extract the 
number of signal S and background B events contained 
in the central cell, by solving the following linear equation 
system: 

Nq = S + B = 76 (7) 
Ni = 0.135- S + 2- B = 51, 
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TABLE VIII: Isoscalar component of the cross section as a function of -Bc.m. for the process e^e —>■ 7^^(1430)/^. Errors are 
statistical only. 



Ec . m . 


(GeV) 


a (nb) 


Ec.m. 


(GeV) 


cr (nb) 


Ec.m.. 


(GeV) 


a (nb) 


Ec.m. 


(GeV) 


a (nb) 


1.92 


-1.96 


0.069 ± 0.036 


2.16 


-2.20 


0.252 ± 0.100 


2.40 


-2.44 


0.014 ± 0.007 


2.76 


-2.84 


0.009 ± 0.009 


1.96 


-2.00 


0.048 ± 0.033 


2.20 


-2.24 


0.115 ± 0.060 


2.44 


-2.48 


0.128 ± 0.066 


2.84 


-2.92 


0.022 ± 0.012 


2.00 


-2.04 


0.083 ± 0.051 


2.24 


-2.28 


0.059 ± 0.037 


2.48 


-2.52 


0.008 ± 0.004 


2.92 


-3.00 


0.065 ± 0.019 


2.04 


-2.08 


0.261 ± 0.088 


2.28 


-2.32 


0.034 ± 0.023 


2.52 


-2.60 


0.030 ± 0.031 


3.00 


-3.08 


0.044 ± 0.025 


2.08 


-2.12 


0.172 ± 0.089 


2.32 


-2.36 


0.045 ± 0.028 


2.60 


-2.68 


0.012 ± 0.012 








2.12 


-2.16 


0.018 ± 0.018 


2.36 


-2.40 


0.063 ± 0.055 


2.68 


-2.76 


0.003 ± 0.003 









TABLE IX; Isovector component of the cross section as a function of Ec.m. for the process e^e — > _?('|(1430)J('. Errors are 
statistical only. 



Ec . m . 


(GeV) 


a (nb) 


Ec.m. 


(GeV) 


a (nb) 


Ec.m. 


(GeV) 


a (nb) 


Ec . m . 


(GeV) 


a (nb) 


1.96 


-2.00 


0.10 ± 0.07 


2.20 


-2.24 


0.21 ± 0.07 


2.44 


-2.48 


0.06 ± 0.03 


2.84 


-2.92 


0.04 ± 0.02 


2.00 


-2.04 


0.23 ± 0.09 


2.24- 


-2.28 


0.42 ± 0.10 


2.48 


-2.52 


0.02 ± 0.01 


2.92 


-3.00 


0.05 ± 0.02 


2.04 


-2.08 


0.08 ± 0.05 


2.28- 


-2.32 


0.32 ± 0.11 


2.52 


-2.60 


0.06 ± 0.03 


3.00 


-3.08 


0.03 ± 0.01 


2.08 


-2.12 


0.19 ± 0.09 


2.32- 


-2.36 


0.15 ± 0.07 


2.60 


-2.68 


0.13 ± 0.03 
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-2.16 


0.33 ± 0.10 
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FIG. 21: The isospin amplitudes (defined as the square root 
of the isospin components divided by the phase space) for the 
KsK'^TT^ channel, obtained by fitting the Dalitz plot in the 
region -Ec.m. > 2.0 GeV, in the case of 7^2(1430)^^ interme- 
diate state. 



where Nq and A^i are the yields in the central cell and the 
sum of yields in adjacent cells (i.e. the cells numbered 
as 2, 4, 6 and 8 in Fig. 23), respectively. The fractions 
of signal events contained in the central and adjacent 
cells have been determined from MC simulation, while 
the factor of 2 multiplying B in the second equation is 
a consequence of the fact that the background is linear 
with respect to both 77 and KK masses. We find 54 
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FIG. 22: (a) K'iK'^ TT^ mass distribution in the J/ip mass 
region. (b-d) Mass distributions of selected events for: 
(b) K*°{892)K'^ {K*°{892) ^ K^tv^)- c) A'*±(892)Ji"F 
(7^*^(892) ^ K'^Tv'^) and (d) K*^{892)K'^ (7^*^(892) 
A"±7r°). The curves represent the fits obtained using a sum 
of two Gaussians and a linear function. 



signal and 22 background events. 

The simulation of specific reactions with a real 



and/or tt*^ in the final state, like e^e 



&7r*'7r'^7, e ' e 



07, e^e 



KKtt^tt^^ (non-resonant KK produc- 
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TABLE X: Branching ratios of the J/ip KKtv decays. The 



value of Fee is taken from Ref. [171 



Final state, / Yield 



S^'*Ff/*(eV) 



L. 



{K^iT^)Ks 94±9 8.85 ±0.85 ±0.50 1.59±0.16±0.09 
{Ks-K^)K'^ 89±9 8.38 ±0.85 ±0.50 1.51±0.16±0.09 
{K^-K°)K^ 155±12 10.96±0.85±0.70 1.97±0.16±0.13 



TABLE XI: Branching ratios of the J/i) ^ K* K decay {vH'^ 
from Ref. [17] and c.c. stands for charge conjugate). 



Bf^- 10» 



Final state, / 

{K*+K-)+c.c. 29.0±1.7±1.3 5.2±0.3±0.2 
(i^'^T^Vc-c. 26.6±2.5±1.5 4.8±0.5±0.3 



• 10^ 
(Ref. [17]) 



5.0±0.4 
4.2±0.4 
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FIG. 23: The 2-dimensional if^ , 77 invariant mass distri- 
bution around the and vr" masses. The central cell identifies 
the signal region. 



tion), and gg, tt production from continuum, shows that 
their contribution is negUgible. We can thus conclude 
that the only sources of background are represented by 
the non-resonant KKtjP^ final state, and by a number 
of ISR processes with a true and a random 77 combi- 
nation. Figure 24(a) shows the signal in the K^K~ 
invariant mass, after the cut around the tt*^ mass is ap- 
plied; viceversa, Fig. 24(b) shows the 7r° signal once a 
is selected. Fitting these distributions with a Gaussian 
and a linear background we estimate the two background 
sources to be 15 ± 2 and 7 ± 3 events, respectively. 

The mass spectra for the data and the estimated back- 
grounds arc shown in Fig. 25. The mass spectrum for the 
KK'K^ background is obtained from the simulated sam- 
ple, while the spectrum for the other ISR backgrounds is 
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FIG. 24: The K (a) and 77 (b) invariant mass distri- 

butions. Signal and side-band regions are shown as vertical 
lines in the plots. 



> 

o 

o 
in 
o 
o 



10 



c 

> 




1.5 



2.5 3 
MC^Ti") (GeV/c^) 



FIG. 25: The invariant mass distribution of the (/)7r" system 
(open histogram) . The hatched and gray histograms show the 
background from non-(^ K~ -n^ and from 77 combinatorial. 



taken directly from the 7r° sidebands. 

The e+e^ ^ (?!)7r° process is simulated considering only 
(j)^ decays, so that the cross section as a function 
of the effective cm. energy is given by 



{Ec.m) 



K+K-) ' 
(8) 

where Ec.-m. is the invariant mass of the ^tt" system; 
4> dN^^o^ is the number of selected 07r° events after back- 
ground subtraction in the interval dEc.m., £(£-0.111.) is the 
detection efficiency obtained from MC simulation (see 
Fig. 26), and B{<f> ^ K^) is the world average value 
of the (p^ decay rate [17]. The energy behavior 

of the cross section is shown in Fig. 27 and reported in 
Table XII. The cross section has a maximum of about 
0.2 nb around 1.6 GeV/c^ 
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TABLE XII: Measurement of the e^e 4>tx cross section as a function of Ec.m.- Errors are statistical only. 



Ec.m. 


(GeV) 


a (nb) 




(GeV) 


cr (nb) 


Eq . m . 


(GeV) 


a (nb) 


Ec . m . 


(GeV) 


a (nb) 


1.20 


-1.30 


0.014 ± 0.016 


1.55 


-1.60 


0.191 ± 0.060 


1.75 


-1.80 


0.015 ± 0.015 


1.95 


-2.10 


0.002 ± 0.004 


1.30 


-1.40 


0.025 ± 0.018 


1.60 


-1.65 


0.087 ± 0.039 


1.80 


-1.85 


0.046 ± 0.027 


2.10 


-2.30 


0.006 ± 0.005 


1.40 


-1.50 


0.033 ± 0.022 


1.65 


-1.70 


0.072 ± 0.035 


1.85 


-1.90 


0.093 ± 0.038 


2.30 


-2.60 


0.002 ± 0.002 


1.50 


-1.55 


0.073 ± 0.038 


1.70 


-1.75 


0.100 ± 0.041 


1.90 


-1.95 


0.089 ± 0.036 
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FIG. 26: Detection and reconstruction efficiency as a func- 
tion of the cm. energy for the 07r° final state. The solid 
line is the result of, a fit to the data with the function 



a 1 



where x = Ec. 
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FIG. 27: The e+e" 



'7r° cross section. 



The efficiency corrections and the systematic uncer- 
tainties which affect the 07r° cross section measurement 
are similar to those of the if+if ^tt" reaction. However, 
the uncertainties on background subtraction are larger, 
and we have to account also for the 1.2% uncertainty on 
the (f> K~^K~ branching fraction. The total system- 
atic error associated to the cross section measurement is 
therefore ±8.7%. 



IX. THE K+K-r] FINAL STATE 

We analyze (firj and K^K~ri (with non-resonant 
K^K~ pairs) final states separately, selecting samples 
with an invariant mass of the K~ system below and 
above 1.045 GcV/c^: and NOT-0 events, respectively. 
In both cases, because of the cleanliness of the 77 signal, 
no additional cuts on photons or extra neutral energy 
are required. The K^K~ and 77 invariant mass distri- 
butions for selected events are shown in Fig. 28. 
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FIG. 28: The K+K' (a) and 77 (b) invariant mass distribu- 
tions for selected events. 



A. e+e" KKri 

Following the same strategy as the analysis described 
in the previous sections, we select a total of 113 candi- 
date events. We estimate 42 ± 42 qq and 15 ± 8 ISR 
background events: 56 ±43 events is then the overall sig- 
nal yield. Figure 29 shows the e+e^ K^K^rj mass 
spectrum from production threshold up to about 4 GeV, 
together with the estimated backgrounds, while the rel- 
ative cross section is depicted in Fig. 30 and the results 
are listed in Table XIII. 
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TABLE XIII: Measurement of the e^e —^ rjK'^K (NOT(/>) cross section as a function of Ec.m.- Errors are statistical only. 



77 f 
-C/c.m.l. 


Kje\ ) 


<7 [XID ) 


-fic.m.l, (jSW ) 


a [no ) 


-C'c.m. 


(GeVJ 


a [no ) 


-fic.m. ( (jCV ) 


G i^nb j 


1.64- 


-1.74 


0.08 ± 0.04 


2.24-2.34 


0.02 ± 0.01 


2.72 


-2.76 


0.02 ± 0.02 


2.96-3.00 


0.03 ± 0.02 


1.74- 


-1.84 


0.01 ± 0.03 


2.34-2.44 


0.01 ± 0.01 


2.76 


-2.80 


0.01 ± 0.01 


3.00-3.04 


0.01 ± 0.01 


1.84- 


-1.94 


0.06 ± 0.02 


2.44-2.54 


0.03 ± 0.01 


2.80 


-2.84 


0.01 ± 0.01 


3.04-3.08 


0.01 ± 0.01 


1.94- 


-2.04 


0.03 ± 0.02 


2.54-2.64 


0.01 ± 0.01 


2.84 


-2.88 


0.01 ± 0.02 


3.08-3.10 


0.27 ± 0.09 


2.04- 


-2.14 


0.03 ± 0.01 


2.64-2.68 


0.01 ± 0.01 


2.88 


-2.92 


0.01 ± 0.01 


3.10-3.12 


0.30 ± 0.09 


2.14- 


-2.24 


0.00 ± 0.01 


2.68-2.72 


0.01 ± 0.02 


2.92 


-2.96 


0.02 ± 0.02 


3.12-3.14 


0.11 ± 0.05 



The uncertainties on background estimation, in par- 
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FIG. 29: The invariant mass distribution of the K~ r\ sys- 
tem (open histogram). The hatched and the gray histograms 
show the qq and ISR backgrounds, respectively. 

ticular the qq contribution, completely dominate the sys- 
tematic errors, which sum up to about 80%. On the 
other hand, the J/i/' region is essentially free from any 
background, as can be seen from the very clean 77 and Jji]} 
signals shown in Fig. 31. We select 21±3 Jjtl) K~^K~r] 
events, from which we measure 



B 



■ T'W = (4.8 ± 0.7 ± 0.3) eV, (9) 



FIG. 30: The e"'"e K^K rj cross-section as a function of 
cm. energy. 
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and extract the first evaluation of the J/ip 
branching fraction: 



Bj/^^K+K-i) 



.7 ± 1.3 ±0.7) 10" 



K+K 



(10) 



FIG. 31: (a) the 77 invariant mass distribution in the 77 region, 
when a J/ip K'^K~rj decay is selected, (b) the final J/^p 
signal (0.513 < M^^ < 0.583 GeV/c^ required). 



B. 



The extremely clean (f> signal makes this channel almost 
background- free. Following the background-subtraction 
procedure as for the (f>Tr'^ final state we get 483 signal 



and 9 background events. The (j> signal, once a 77 is se- 
lected, and the 77 signal, after the (j) selection, arc shown 
in Fig. 32. 

We simulate several specific reactions with real cf> 
and/or j] in the final state, namely different ISR pro- 
cesses and qq, tt production, finding no significant con- 
tribution. The mass spectrum for signal events is shown 
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in Fig. 33, together with the subtracted distribution for 
the background events obtained from the side bands. 
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FIG. 32: The K~ (a) and 77 (b) invariant mass distri- 
butions. Signal and side-band regions are shown as vertical 
lines in the plots. 
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FIG. 33: The invariant mass distribution of the (j>ri system. 
The shaded histogram refers to the small subtracted back- 
ground component. 

The efficiency corrections and the systematic uncer- 
tainties that affect the (pr] cross section measurement are 
very similar to those of the K^K~tt'^ reaction (see Ta- 
ble IV). However, the background level is so low for the 
(j)rj final state that the associated errors are negligible 
and do not affect the background subtraction procedures. 
The systematic uncertainty associated with rj reconstruc- 
tion is equivalent to that for photon reconstruction, and 
the uncertainty on the (j) K^K~ decay rate is 1.2%. 
The total systematic error is thus ±5.3%, while the effi- 
ciency is overestimated by 4.6%. The detection efficiency 
shown in Fig. 34 as a function of cm. energy, is obtained 
from MC and accounts also for 77 — > 77 decay rate. The 
corresponding values for the e^e^ — s- (jyrj cross section. 



FIG. 34: Detection and reconstruction efficiency as a func- 
tion of the cm. energy for the (jyrj final state. The solid 
line is the result of, a fit to the data with the function 
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FIG. 35: e^e~ (j>ri cross section distribution as a function 
of cm. energy. 



obtained with an expression analogous to Eq. 8, are re- 
ported in Table XIV and in Fig. 35. 

The cross section peaks at about 3.0 nb around 
1.7 GeV and also hints at a possible resonant behavior 
around 2.1 GeV. Above the J/ip mass the cross section 
is very small, as shown in the inset of Fig. 35. 
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TABLE XIV: Measurement of the e^e — » (^77 cross section as a function of -Ec.m.- Errors are statistical only. 
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X. CROSS SECTION FITS 



The cross section parametrization used to fit the data 



leading to the final state / 



/, at s 



Ay- 



is)- 



E 

R 



, (11) 



where Vfis) is the phase space of the / final state, 
/l" '^ (s) describes the non-resonant background, mainly 
due to the tails of resonances below threshold, and the 
sum runs over all the vector resonances, with mass M^, 
width T fi and relative phase r, assumed to contribute 
to the cross section. All of the final states analyzed eon- 
tain a vector and a pseudoscalar meson. The phase space 
for / = VP has the form: 

P.H^)=[(- + <-^^)^-^<-1\ (12) 



A. Fitting KKn and (^77 cross sections 

The statistics arc large enough to attempt a complete 
description of the A' A'* (892) channel. The K^^K^ti^ 
Dalitz plot yields not only the isoscalar ctq and isovec- 
tor (Ji components, but also their relative phase A0 
(Figs. 16, 17 and 18). In addition, the e+e" 
if±if*T(892) 

cross section, obtained by integrating the 
K'^K~'k'^ Dalitz plot multiplied by the isospin factor of 
three, and taking into account corrections due to the in- 
terference relevant for K^n^ in the K*{^%2) region, pro- 
vides a further constraint on ctqi o"! A0. Indeed, 
both the cross sections for e+e^ K^K*^{^%2) and 



e+e- ^ i^OA' (892)+c.c. can be defined in terms of 
opposite combinations of the same isospin components: 



K'>K (892) +C.C 



(Tie 



iA<t> I 



(13) 



(To + <Ti + 2 Y^owi COS A0 



0'K±if=F(892) 



(Tie 



= aQ + ai — 2-y/(To(7i cos A<j). 



We note that the square root of the cross sections uo 
and (Ti, having roughly the same phase space (the mass 
difference between neutral and charged AT* (892) can be 
neglected here and in the following), are proportional to 
the corresponding amplitudes. The isoscalar A' A'* (892) 
cross section co and cr^^ appear to be dominated by the 
same structure, a broad peak with M « 1.7 GeV/c^. 
This behavior is theoretically expected because the two 
final states have the same quantum numbers I'^{J^'~^) = 
0~(1 ) and similar strangeness content. In principle, 
quantum number conservation should allow w-like and 
0-like contributions for both cross sections. However, in 
the case of the (prj final state, the OZI rule [23] strongly 
suppresses any intermediate w-recurrence. On the other 
hand, for the channel A' AT* (892), where the OZI rule does 
not apply, there are also suppression factors. Suppression 
of w-like versus (/)-like coupling could account for as much 
as a factor of four [17, 24]; an additional price has to 
be paid if a strange quark pair has to be created from 
the vacuum. Hence, in the following we will consider all 
the structures contributing to the isoscalar amplitudes 
as 0-like resonances. In particular, the peak with M w 
1.7 GeV/c^ matches the first 0-like excited state, so we 
will assume it to be simply the 0'. 

The first result that can be drawn comparing the cross 
section data on ctq and tr^^ (Figs. 16 and 35) is an esti- 
mate for the ratio between the decay fractions of the 
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into KK* (892) and 0?7. Indeed, at the (j)' peak: 



E 



'^KK* {892) 



(14) 



where is the branching ratio for the decay R ^ f . 
From Eq. 14 one can infer that both the A'/v*(892) and 
the 077 channel represent an important fraction of the 0' 
total width. 

In order to make optimal use of the experimental data, 
we perform a global fit exploiting all available informa- 
tion from the KKir final states (fom' sets of data: ctqj <^1i 
Acj) and cr^±^.T(8g2)) as well as the 4>r] cross section. 

We have considered also a recently published BaBAR 
measurement on the (j)r] cross section [25] in the e+e^ ^ 



K+K- 



'TT^ channel. These data, even though with a 



lower statistics, are in agreement with our measurement. 

We have performed two fits, in the first case we 
used only our samples (five sets of data: tro, cri, A0, 
<^K±K*T (892) and cr^^), in the second case we included 
the (pr] cross section from Ref. [25] . 

We parametrize the KKtt isoscalar and the isovector 
cross sections using Eq. 11 with one resonance each and 
a non-resonant background, assumed to be the (f) tail. 
The isoscalar cross section shows a clearly resonant be- 
havior, the shape of the sub-dominant isovector cross sec- 
tion instead appears compatible with a pure phase space. 
However, using the information of the relative phase, Acjj 
(Fig. 18), it can be shown that such a quantity is not 
consistent with a real isovector amplitude, as expected 
in case of a non-resonant behavior. If this were indeed 
the case, the relative phase should be as shown by the 
light gray band in Fig. 36. This is clearly not compatible 
with the data (solid circles). The isovector cross section 
tTi , then, has to be described in terms of a p-like excited 
state, the p' in the following. 

Finally, the (j)?] cross section is parametrized using two 
resonances: in addition to the dominant 0' just above 
threshold, we include a small contribution, probably due 
to an excited ss or a more exotic state (tetraquark. hy- 
brid) with a mass around 2.13 GeV/c^, the 0", that has 
been observed decaying into 0/o(98O) [26]. 

A good parametrization of the non-resonant back- 
ground, the function Ay-{s) (/ = KK* (892), cjjr]), can 

be obtained using a 1/s scaling (yl^ '^-(s) A'j-'' /s). As 
a result, the shape of the non-resonant part, Eq. 11, is 
given directly by the phase space, and we may define the 
non-resonant cross section as: 



12nVf{s) [Ay-fs] 



(15) 



We stress that, as shown in Table XV, this cross section 
represents only a small fraction of the resonant one. 

We describe wide resonances using energy-dependent 
widths, Tfi(s) (see Eq. 11), parametrized as: 



rfl(s) 



(16) 



where Tn is the constant total width, and Vk{s) and 
are the phase space and the branching fraction for 
the transition of R into the final state k respectively 
{J2k^k = ca^*^ narrow peaks we use fixed 

widths. 
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FIG. 36: Phase difference between isovector and isoscalar 
KK* (892) amplitudes: data and fits superimposed (dark-gray 
and hatched bands). In particular the dark-gray band is ob- 
tained by using only our (pi] cross section, while the hatched 
band including also the data from Ref. [25]. The light-gray 
band represents the relative phase expected in case of non- 
resonant isovector component. 



As already mentioned, the most important resonance 
in these processes is the (j)': it represents the main con- 
tribution for both ctq and cr^,,. It follows that the total 
width T^r{s) will account for both these VP-like final 
states: the isoscalar KK* (892) and the (prj. Other final 
states with small branching ratios are parametrized by 
adding a constant term to the width 



r0'(5) 



r 



KK- (892) 



■Bl 



VKKH892){Ml,r^^'('''^ 



(17) 



'^KK'{892) '^<j>ri 



The <f>" , with a constant width, is included to model the 
077 cross section. 

The main contribution to the cross section cri comes 
from the p' resonance. We assume the p' to decay mainly 
in four pions, and again include a constant term to ac- 
count for other channels: 



rp'(.s) = r. 



, (18) 
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FIG. 37: Isoscalar (a) and isovector (b) components of the 
KsK^TT^ cross section. The sohd points are the data and 
the bands represent the fits. The hatched band is obtained 
including the (pi] cross section from Ref. [25]. 



where, for Vi-„{s) we use: 



(19) 



with n = 3/2. Note that the fit is not sensitive to this 
power, as higher values (e.g. n = 5/2) give similar pa- 
rameters. 

Our global fit results are summarized in Table XV and 
shown in Figs. 36-39, where the errors of data points in- 
clude both statistical and systematic uncertainties, and 
the bands represent the convolution of the one-sigma sta- 
tistical errors of the fit. The solid and the hatched fill 
refer to the two considered cases. In particular. Figs. 37 
and 36 show data and fits for the isospin components and 
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Ec...(GeV) 

FIG. 38: Data on the annihilation cross section cr^±_ff*T{892) i 
see Eq. 13, and descriptions (solid and hatched bands) in 
terms of the parameterization of Eq. 11, with parameters re- 
ported in Table XV. The hatched band is obtained including 
the (^77 cross section from Ref. [25]. 




FIG. 39: Data on the annihilation cross section a^r) (the solid 
circles are from this analysis, while the empty squares are from 
Ref. [25]) and descriptions (solid and hatched bands) in terms 
of the parameterization of Eq. 11, with parameters reported 
in Table XV. The hatched band is obtained including the 077 
cross section from Ref. [25]. 



the relative phase from the ifgi^^Tr^ Dalitz plot anal- 
ysis, Fig. 38 shows the cross section cr^± if (892): and 
Fig. 39 shows cr^^, as obtained in this analysis and in 
Ref. [25]. 



25 



TABLE XV: Parameters for resonances (f)' and p' obtained by fitting simultaneously the isoscalar and isovector component of 
the KK*{9,^2) + c.c. cross section, their relative phase, the K'^ K~-k^ and the (^77 cross sections (only our samples on the left, 
including also data from Ref. [25] on the right). Real backgrounds are defined through Eq. 15. All phases r (see Eq. 11) are 
set to zero. 



= 1.28 



R with 1 = 



4>' 



r e&'^KK" (892) 

r«Sj,(eV) 

1 R« 

'-'if if * (892) 

MH(MeV) 
rH(MeV) 



(eV) 



^bkg 



(892)(A^|')(nb) 



367±47 

154±32 1.7±0.8 
0.33±0.14 

1709±19 2127±24 
325±68 60±50 
0.8±0.3 

(4.7±1.4) X 10"^ 



R with 1 = 1 



P 



ee'^KK* (892) 

1 - sf. 

MH(MeV) 
rfl(MeV) 



(eV) 



.(892)(^4)(nb) 



129±15 
0.36±0.22 
1508±19 
418±26 
0.9±0.2 



= 1.20 



369 ±53 

138±33 1.7±0.7 
0.36±0.15 

1709±20 2125±22 

322±77 61±50 
0.9±0.3 

(4.7±1.8) X 10"-' 



127±15 
0.35±0.23 
1505±19 
418±25 
0.9±0.2 



All these data sets support the hypothesis that the 
main contribution in the isoscalar amplitudes is the res- 
onance (j)' , while the isovector component is almost fully 
described by a broad excited p state. However, the nor- 

2 

malized ^ values, 1.28 in the first case and 1.20 in the 
second case when the data on the 0?/ cross section from 
Ref. [25] are used, do not exclude the presence of other 
possible structures but, on the other hand, the accuracy 
of the data does not allow any clear identification. 

As mentioned above, the mass and width of the are 
compatible with the </)(1680) [17]. The parameters ob- 
tained for the (/)", the second isoscalar excited state, arc 
consistent with those reported in Ref. [26]. A clear evi- 
dence for this decay channel cannot be established with 
the available data sample, however, the probability to be 
a statistical fluctuation is 6x10"'^ (25 events under the 
peak with an expectation of 11 evaluated from the mass 
spectrum in Fig. 35). The p' from the isovector compo- 
nent of the KK* {892) cross section has a mass value in 
agreement with Ref. [17], although this agreement is not 
particularly striking given the width of this resonance. 

It is worth stressing that at i?c.m. above 1.8 GeV, con- 
tributions coming from any higher mass K* certainly af- 
fect the fit in the KK* (892) channel. At low Ec.n,. the fit 
is very sensitive to the if* (892) width, when the Ktt in- 
variant mass is close to the K*{892) threshold. However, 
analyticity relates if i^* (892) and (j)r] amplitudes to the 
i4:*(892) and <j) radiative widths (that is A' A'* (892) and 
(j)ri at s = 0) via dispersion relations. Exploiting these 
relations might lead to an improved determination of the 



threshold cross section. 

The systematic effects due to the choice of fit model 
were evaluated by repeating the fits using different val- 
ues for the relative phases '^r, and by choosing different 
parameterizations for the energy-dependent widths, as 
appropriate. 



B. Fitting the 077° cross section 

The (pir'^ channel is a pure isospin one final state, 
hence, using Eq. 11 to describe the corresponding cross 
section, only isovector contributions must be included. 
This channel has been selected as a possible candidate 
for multiquark intermediate states: contributions are ex- 
pected from p recurrences, even though they are OZI- 
suppressed. We label any resonance in this channel as a 
/9-like excited state. 

Two cases have been studied in detail: 

• only one p-like excited state, a p" (to distinguish it 
from the p' of the previous case); 

• two p recurrences, i.e., in addition to p", a nar- 
row resonance to better represent the data around 
Ec.m. « 1900 MeV (Fig. 40), labelled in the follow- 
ing as /9(1900). 

The p(1900) width is assumed constant, while for the p" 
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TABLE XVI: Parameters obtained for the <f>Tv'^ cross section. 
First column: with the only p" resonance, second and third 
columns: including also the /9(1900). The normalized a-nd 
corresponding C.L. are reported in each case. 



(^,C.L.) 


(#^.0.31) 


( 7.37 
V 16-8 ' 


0.50) 


R 


P 


P 


p(1900) 


r«Z?J^o(eV) 


4.4±1.0 


3.5±0.9 


2.0±0.6 


(1 - ) 


0.67±0.43 


0.44±0.49 




MH(MeV) 


1593±32 


1570±36 


1909±17 


rfl(MeV) 


203±97 


144±75 


48±17 


*fl(rad) 








TV 


%to(-'W'p(i9oo) 


)(nb) (0.4±0.2)xl0"-' 


(0.5±1.5)xl0"=' 



(GeV) 



XI. 



SUMMARY AND CONCLUSIONS 



FIG. 40: Data on the annihilation cross section cr^^o and 
descriptions in terms of the parameterization of Eq. 11 with 
the only p" , hatched band, and including also the p(1900), 
gray band. 



we use the energy-dependent width 



K>P 



4-K 



^ "47r 



where 7^4^ (s) is the four-pion phase space defined in 
Eq. 19, and is the branching fraction for p" Air. 

The results of the fit with only the p" are reported in 
the first column of the Table XVI. 

In the second case, where the p(1900) is also included, 
a two step procedure is used to determine the phases of 



the two quasi-real amplitudes. We find (^p" , ^ 



p(1900) 



(OjTt) to be the best combination. With these phases we 
get the results reported in the second and third columns 
of the Table XVI. 

The results of the two fits are shown in Fig. 40, super- 
imposed on the cross section data. 

A slightly better is obtained by adding this extra 
resonance. We can not however exclude that the observed 
accumulation of events at i?c.m. ~ 1.9 GeV is produced 
by a statistical fiuctuation. In fact, we observe 18 events 
to be compared to an expectation of 8 events, taking into 
account both the background and the tail of the p", and 
this translates in a Poisson probability of 2 x 10^'^. 

Mass, width, and quantum numbers [I^{J^'"') = 
1+(1~~)] obtained for the p(1900) are compatible with 
those of the so-called "dip" observed in other channels, 
primarily multi-pion final states [27] . 



We have studied the K^K^tt"^, K+K-t:°, 07rO, 
K~^K~T], and (prj final states, produced in BaBAR via ISR. 

We have measured the production cross section from 
threshold up to «4.6 GeV for all of these processes to an 
unprecedented accuracy. 

By studying the asymmetric K^K^n^ Dalitz plot, 
we have obtained the moduli and relative phase of the 
isospin components for the KK* (892) cross section. The 
knowledge of the isoscalar and isovector cross sections, 
(To and (71, allows a simple description in terms of vec- 
tor meson resonances, where p- and 0-like mesons do not 
mix. In addition, the relative phase gives unique infor- 
mation on the sub-dominant isovector component, which, 
although compatible with a pure phase space behavior, 
reveals a resonant structure. 

Two global fits, that benefit from five interconnected 
sources of information (cto, ci, cr^,;, (Jk+k-ttO ^-nd A0), 
have been performed. In a first case we have used only the 
data samples coming from this analysis, while in the sec- 
ond case recent BaBAR data on cr^,, in the A'+if ^tt+tt^tt" 
channel have been included. 

Both (prj and the isoscalar A' A'* (892) cross sections 
have been parametrized with the same dominant reso- 
nance 0'; the suitable combination of ctq, (Ti and A<^, 
obtained from the A'^AT^tt^ channel, is used to describe 
the Ar+Ar~7r° cross section. The mass and width for 
the 0', the main resonant contribution in the isoscalar 
channels, are compatible with those of the (/)(1680) me- 
son reported in Ref. [17]. Concerning the (p" resonance, 
included only in the (f)?] channel, we find parameters com- 
patible with those of the first observation in the 0/o(98O) 
final state [26] , confirming the quantum numbers of this 
resonance, I^{J^^) = 0"(1 ). 

The isovector component of the AT AT* (892) cross sec- 
tion is described by one broad resonance, p', whose mass 
is compatible with the p(1450), the first p recurrence. 
The slight inconsistency for the width is probably due to 
the fact that a single structure is not sufficient to describe 
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TABLE XVII: Summary of parameters of resonances in the studied final states, obtained including also the data on the (prj 
cross section from Ref. [25]. The parameters for the p" are taken from the fit to the (^7r° cross section with two resonances. 



Isospin 


R 


rfeSjJjf.(g92)(eV) 


rf,i3^,(eV) 


Mfl(MeV) 


rfl(MeV) 





<P' 


369±53±1 


f38±33±28 


f709±20±43 


322±77±160 






1.7±0.7±1.3 


2125±22±f0 


61±50±f3 


isospin 


R 


^ee^KK* (892) ( ^) 


r«S«o(eV) 


AfH(MeV) 


rfl(MeV) 




1 

P 


127±15±6 




1505±19±7 


418±25±4 


1 


P" 




3.5±0.9±0.3 


1570±36±62 


144±75±43 




p(1900) 




2.0±0.6±0.4 


f909±f7±25 


48±17±2 



this cross section and that the fit tends to broaden the p' 
to mimic a more complex structure. The present statis- 
tics does not allow the identification of multiple struc- 
tures. 

The (jji:'^ cross section has been fit with p recurrences, 
whose coupling with the final state is, however, OZI- 
suppressed. Two possible descriptions have been con- 
sidered, with one and two resonances. Their parameters 
are reported in Table XVI. 

The resonance labeled p", might be the so called 
C(1480), observed in -K^p (pw'^n charge-exchange re- 
actions [28], in the same (/)7r° final state. However, a 
firm conclusion cannot be drawn at the moment; an OZI- 
violating decay of the meson p(1700) [17] cannot be ex- 
cluded. The second structure, p(1900), is compatible 
with the "dip" already observed in other experiments, 
predominantly in multi-pion final states [27] . 

In Table XVII we summarize the whole set of parame- 
ters obtained from the fit to the cross sections, with the 
inclusion of model related systematic errors. 

As already pointed out, the normalized x^'s obtained 
from this global-fit procedure, are still compatible with 
other possible minor contributions that, however, can- 
not be precisely identified because of the accuracy of the 
data. 

In the (2 — 3) GcV cm. energy region the K'^K^tt^ 
final state has been analyzed, adding the K2{14:30)K to 
the K*{892)K intermediate state. The strong asymme- 
try between neutral and charged channels of the Dalitz 
plot (Af^^^i vs. M^±^:p) can be explained in terms of 
constructive and destructive interference between differ- 
ent isospin components. This asymmetry might be con- 
nected to a similar effect observed in the radiative decay 
rates of the neutral and charged iir|(1430) [17]. 

Very clean J/ip signals have been observed in most of 
the studied final states, allowing the measurement of the 
corresponding branching fractions. 



We have no evidence of the y(4260) decay- 
ing into KgK^n^ and K'^K^tt'^. We determine 

rl^''''^BlfJ'!l < 0.5 eV and vl^''''^ bI^k^o < 0.6 
cV at 90% C.L.. We observe also (inset of Fig. 6) a 
3.5 (T fluctuation in the e+e~ ^K'^K^n^ cross section 
at i^c.m. ~ 4.20 GcV. 
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